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Abstract 

We investigate theoretically a Fano interferometer composed by STM and AFM tips close to a Kitaev dimer 
of superconducting adatoms, in which the adatom placed under the AFM tip, encloses a pair of Majorana 
fermions (MFs). For the binding energy A of the Cooper pair delocalized into the adatoms under the tips 
coincident with the tunneling amplitude t between them, namely A = t, we find that only one MF beneath 
the AFM tip hybridizes with the adatom coupled to the STM tips. As a result, a gate invariance feature 
emerges: the Fano profile of the transmittance rises as an invariant quantity depending upon the STM tips 
Fermi energy, due to the symmetric swap in the gate potential of the AFM tip. 
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1. Introduction 

The physicist Ettore Majorana proposed in the 
field of high-energy Physics the existence of pecu¬ 
liar fermions that constitute their own antiparticles. 
In the context of condensed matter Physics, these 
fermions are Majorana quasiparticles 0,0. From 
the quantum computing perspective, two Majorana 
fermions (MFs) can compose a regular fermion act¬ 
ing as a protected qubit, which is indeed decoupled 
from the host environment and free of the decohcr- 
ence effect. As a result, the quest for setups sup¬ 
porting MFs has attracted broad interest from the 
communities of researchers in the fields of quantum 
information and transport M, since the qubit 
made by the coupled MFs appear only in the topo¬ 
logical phase. Noteworthy, the Kitaev chain within 
such a phase [6j is considered the most promising 
candidate to this end as the aftermath of the emerg¬ 
ing p -wave and spinless superconductivity. Indeed, 
in Kitaev’s setup, MFs appear as zero-energy modes 
attached to the edges of the chain. 
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Experimentally, p -wave superconductivity is fea¬ 
sible due to the proximity effect by the employ¬ 
ment of an s-wave superconductor close to a semi¬ 
conducting nanowire characterized by a spin-orbit 
interaction under an external magnetic field @- 
0- Particularly in the case of transport through 
quantum dots (QDs) coup led to a MF fllFf23| . a 
zero-bias peak (ZBP) [l5|, [l6| in the conductance 
is expected to be observed. It is worth mention¬ 
ing that the ZBP has been detected in conductance 
measurements through a nanowire of indium an- 
timonide linked to gold and niobium titanium ni¬ 
tride electrodes [24]. Analogously, a ZBP has also 
been verified in the superconducting system of alu¬ 
minium next to a nanowire of indium arsenide [251. 
However, the ZBP signature may also have another 
physical origin, for instance the Kondo effect [26 
1301 1 ■ thus turning the experimental ZBP detection 
inconclusive within a MF perspective. Moreover, 
recently an alternative way for the achievement of 
the topological Kitaev chain has been the employ¬ 
ment of magnetic chains on top of superconductors 
IrD-Ffill. Particularly in Ref. 341, the ZBP observed 
exhibits a subtle amplitude of the order 10 _4 (2 Gq), 
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Figure 1: (Color online) (a) Setup composed by STM and 
AFM tips in the presence of superconducting adatoms hosted 
by a conventional superconductor with strong spin-orbit cou¬ 
pling. (b) The adatom 2 is made by a pair of coupled Ma- 
jorana fermions (MFs A and B represented by half-spheres) 
where the level £2 is induced by the AFM tip and plays the 
role of the connection between the MFs. It can assume neg¬ 
ative values when it stays below the MF zero mode, due to 
the tuning performed by the gate potential of the AFM tip. 
(c) In the case of i 7 ^ A and £2 7 ^ 0, both MFs hybridize 
with the adatom 1. (d) For the matching of the Cooper 

pair binding energy with the normal tunneling between the 
adatoms 1 and 2 (t = A) and £2 7 ^ 0, the MF A becomes 
the unique Majorana connected to the adatom 1. (e) Here 
we see the isolated MF A coupled to the adatom 1 when 
t = A and £2 = 0, which results in the standard ZBP in the 
transmittance and in the isolated MF B at one edge of the 
Kitaev dimer. 


which is a signal extremely weak due to thermal 
broadening, where Go = e 2 /h is the quantum of 
conductance. Thus in the current context, novel 
approaches in the pursuit of MFs become necessary. 
In this scenario, we highlight the proposal found in 
Ref. 37|, which is a pioneering example concerning 
this issue: a hybrid spin-microcantilever system for 
optical detection of MFs on the edges of a Kitaev 
chain. 


In this work we propose a new route for detection 
of MFs signatures. To that end, we consider the 
setup outlined in Fig. |T] where an s-wave supercon¬ 
ductor with strong spin-orbit coupling hosts a set 
of magnetic adatoms, in analogy to the experimen¬ 
tal apparatus developed in Ref. [34J that describes a 
superconducting surface of lead (Pb) well known by 
its strong spin-orbit coupling, thus allowing p -wave 
superconductivity on this chain due to the prox¬ 
imity effect. It is worth mentioning that the pro¬ 
posal of Ref. [34| is distinct from the semiconducting 
spin-orbit wire reported in Ref. 24]. In the current 
work, we consider for a sake of simplicity, just a 
pair of superconducting adatoms (Kitaev dimer), 
in particular placed nearby STM (Scanning Tun¬ 
neling Microscope) and AFM (Atomic Force Mi¬ 
croscope) tips. Additionally, we assume two highly 
spin-polarized STM tips in order to freeze the spin 
degree of freedom, thus avoiding the emergence of 
the Kondo effect within the adatom 1. 

In regard to the possibility of experimental re¬ 
alization of the setup shown in Fig.[T| we stress 
that multitip STM experiments can be possible, see 
e.g. Refs. 38|-4(|. In the case of the AFM tip, 
which is fixed on top of an adatom, it operates 
similarly to the Scanning Gate Microscopy (SGM) 
technique [41], wherein a charged tip allows that its 
gate potential tunes the energy levels of the adatom 
probed. This approach is widely employed in sys¬ 
tems with quantum point contacts (QPCs) formed 
in two dimensional electron gases [42 l If.'II , which we 
invoke here to our proposal. Thereby despite the 
challenging of applying the SGM technique in the 
scenario of Fig.[Q we trust that in the near future 
such a procedure can be implemented. 

In the frame of the setup here proposed, a device 
based on Fano interference [H Mm Si becomes 
an alternative method for detection of MFs. As 
a matter of fact, the Fano effect is found in sys¬ 
tems where distinct tunneling paths compete for 
the electron transport. Noteworthy, such a phe¬ 
nomenon can be perceived experimentally by the 
STM just by measuring the conductance Ha-Sol- 
In the setup of Fig[T]here proposed, the Fano effect 
arises from the adatom 1 coupled simultaneously to 
the pair of STM tips, once the electrons can tun¬ 
nel between such tips or (and) directly through this 
adatom as expected for realistic experimental con¬ 
ditions. Thus depending on the Fano lineshape of 
the conductance, we can identify which path rules 
the quantum transport or if they compete on the 
same footing. Later on, we will see that the Fano 
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parameter of interference qb then dictates the Fano 
profile making explicit the dominant path in the 
system. 

Here we show that the two coupled MFs A and B 
within the adatom under the AFM tip (see FigQ} 
lead to a gate invariance feature in the transport 
experiment. Particularly, when the AFM tip swaps 
symmetrically its gate potential around the MF 
zero mode, it reveals a universality signature in the 
Fano profile of the transmittance. More specifically, 
we find then distinct situations in which the trans¬ 
mittance shares the same profile as a function of 
the Fermi level of the STM tips, in particular when 
the binding energy of the Cooper pair split into the 
adatoms under the STM and AFM tips, is in res¬ 
onance with the tunneling term between them. In 
this case, we show that when the MF A couples 
solely to the adatom 1, the gate invariance char¬ 
acterized by the aforementioned universality is re¬ 
vealed by connecting MFs A and B. 

This paper is organized as follows: in Sec. [2] 
we develop the theoretical model for the system 
sketched in Fig.[T]by deriving the expression for the 
transmittance through such a device (see also the 
Appendix) together with the system Green’s func¬ 
tions. The results are present in Sec.|3l and in Sec. [4] 
we summarize our concluding remarks. 

2. The model 

2.1. System Hamiltonian 

Here we consider the system outlined in Fig.QJa) 
for a pair of highly spin-polarized STM tips con¬ 
nected to the hybrid setup composed by an AFM 
tip and the “host+superconducting adatoms” set. 
Thus such a spinless model reads: 

H = ’^^iakC^f.Cak + £ jdjdj 
otk j 

a.k 

+ (V 12 ^ ' c Ifc c 2 g T tdid\ + Acrfjdi + H.c.), (1) 

kq 

where the electrons in the STM tip a = 1,2 (STM 
tip 1 and 2, respectively) are described by the op¬ 
erator cL (c a fc) for the creation (annihilation) of 
an electron in a quantum state labeled by the wave 
number k and energy s a k = £/. — p, a , with m„ as 
the chemical potential. Here we adopt the gauge 
M 2 = A/i and mi = — Am, with M2 — Mi = 2 Am = eip 


as the bias between the tips, being e > 0 the elec¬ 
tron charge and ip the bias-voltage. Consequently, 
the transmittance through the setup is a function 
of the Fermi energy e = mi = M 2 of the STM tips, 
where the point e = 0 corresponds to the MF zero 
mode. For the adatoms, dj ( dj ) creates (annihi¬ 
lates) an electron in the state £j, with j = 1,2. V 
stands for the hybridizations between the adatom 1 
and the STM tips. These couplings are considered 
the same to avoid Andreev currents in the system 


111 . [12] (see also the Appendix). V 12 is the STM 
tip 1-tip 2 coupling, which will ensure the precise 
renormalization of the energy level in the adatom 1 
as we will see. 


The s-wave superconductor with strong spin- 
orbit coupling enters into the model via the param¬ 
eters t and A, which respectively yield the hopping 
term between the adatoms 1 and 2, and also the 
binding energy of a delocalized Cooper pair split 
into such adatoms: as a net effect of the arrange¬ 
ment of Fig. Oja) considered, p -wave superconduc¬ 
tivity is induced on the pair of adatoms (Kitaev 
dimer) close to the tips, similarly to that reported 
in Ref. [2] for the case of a long Kitaev chain. As 
the adatom 1 is placed within the region separat¬ 
ing the STM tips and the superconducting surface, 
we can safely neglect the tunneling from the former 
into the latter as an outcome of the setup adopted. 
We stress that the role of the AFM tip is indeed 
to gate overlap the MFs enclosed by the adatom 2. 
This clue between the MFs found within such an 
adatom will be clarified later on as the own energy 
level of this adatom, which consists in a degree of 
freedom completely tunable by the AFM tip. Thus 
for a sake of simplicity, we clamped these features as 
the most relevant from the proximity problem be¬ 
tween the s-wave surface and the adatoms. Alterna¬ 
tively to this end, an approach concerning ab-initio 
description on the proximity issue can be helpful, 
but it does not belong to the current scope of the 
paper. 


Here we express the adatom 2 in the Majo- 
rana basis by following the transformation g ?2 = 
^('F.a + i^b) and d\ = ^=(\E 'a - ^b), in which 

T] = \Eq (l = A, B) characterizes a MF operator, 
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G = G 0 J ds(-^fjne), 


(3) 


thus yielding 

"Hmfs = £ 2 d\d 2 + {td\d\ + Ad\d\ + H.c.) 

= i£ 2 ^A^B + ^ (di - d\)^ A 

+ i^-j=-^-(di + d\)^B + ( 2 ) 

where the first term in the second line makes ex¬ 
plicit the MFs within the adatom 2 and the others 
stand for the connections between MFs A and B 
with the adatom 1 (see FigsQJc), (d) and (e)). It is 
worth noticing that the standard Majorana ZBP oc¬ 
curs for t = A, £i ^ 0 and £ 2 = 0, with the conduc¬ 
tance G = 0.5Go due to the MF 71 = i(d\ — d\)/y/2 
existing in the adatom 1 given by the second term 
of Eq.© (— 2 it 7 i'F J 4 ), which has a component on 
the adatom 2, where a completely localized MF is 
observed according to M. Leijnse and K. Flensberg 
in Ref. 0. 

In this scenario, we explore the model regimes 
t = A and t ^ A combined with £ 2 > 0 and £ 2 < 0 
varied symmetrically with respect to the MF zero 
mode. We stress that for a fixed A, the tuning of 
the parameter t can be performed experimentally 
just by changing the distance between the adatoms 
1 and 2, while the AFM tip controls the level po¬ 
sition £2 of the adatom 2 , since such a tip oper¬ 
ates as a gate potential. Particularly for t = A, 
we can verify in Eq. © that the adatom 1 decou¬ 
ples from the MF B (here sketched by the most 
right half-sphere in Fig.QJd)) thus allowing an ex¬ 
clusively connection between this adatom and the 
MF A (the most left half-sphere) as the Hamilto¬ 
nian Hmfs = £ 2 {i^A^B + \) + v / 2A(di - d\)^ A 
points out. For £ 2 ^ 0, the Fano profile of the trans¬ 
mittance becomes invariant under the symmetric 
change of £ 2 with respect to the MF zero mode, 
in contrast with the case f / A where this uni¬ 
versality is prevented. Thereafter, it gives rise to 
distinct transmittance profiles strongly dependent 
on the sign of £2 away from the point t = A. 

Thereby we recognize this invariance as a Majo¬ 
rana property of the system when it is driven to the 
point t — A with symmetric swap of £ 2 . 

2.2. Conductance 

In what follows we derive the Landauer-Biittiker 
formula for the zero-bias conductance G (46]. Such 
a quantity is a function of the transmittance T (e) 


as follows: 


where fp stands for the Fermi-Dirac distribution. 

We begin with the transformations c 2 k = 
^75 ( c e k+c 0 k ) and Ci fe = g(c efe -c ofe ) on the Hamil¬ 
tonian of Eq. ©, which starts to depend on the 
even and odd conduction operators c e k and c 0 k , re¬ 
spectively. These definitions allow us to express 
Eq. © as "H = 7i e + Ho + Tdtun = ’% ¥ >=o + Tdtun, 
where 

T~Le — ^ v - A:Cf:A: A £\d-^d\ -f- Fl2 ^ ' C e f~C e p 

k kp 

+ V2V^l k d 1 +U.c.) + n M F S (4) 

k 

represents the Hamiltonian part of the system cou¬ 
pled to the adatoms via an effective hybridization 
V2V, while 

n ° = Yl £ k c lk c ok - Vj 2 C lk C op (5) 

k kp 

is the decoupled one. However, they are connected 
to each other by the tunneling Hamiltonian Jitun = 
~A-t 1 J2k( C ek C ° k + C lfc C efe)- 

As in the zero-bias regime Ap —> 0, due to p —> 0, 
"Htun is a perturbative term and the linear response 
theory (see the Appendix) ensures that 

T (e) = ( 2 ttVi 2 ) 2 p e (£)p 0 (£), ( 6 ) 

where p e {c) = — is the local density of 

states (LDOS) for the Hamiltonian of Eq. © and 

= - jd (r) Tr{£> e [^ e (t) , T* (0)] + } (7) 

gives the retarded Green’s function in the time 
domain r, where 9(r) is the Heaviside step func¬ 
tion, g e is the density-matrix for Eq. ©, T e = 
/ e + {■KYp 0 ) l Aqd 1 is a field operator, with f e = 
J2 p c ep , the Anderson parameter T = 27rF 2 po, with 
po as the density of states for the STM tips and 
q = (ttpoIT 1 / 2 (*g) • 

To calculate Eq. © in the energy domain e, 
we should employ the equation-of-motion (EOM) 
method [46j summarized as follows 

(e + tO+J&is = [A, B'] + + g [WB ( 8 ) 
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for the retarded Green’s function Gab = 
f drGABE^ e+M ) r , which is the time Fourier trans¬ 
form of GaBi with A and B as fermionic operators 
belonging to the Hamiltonian Hi. By considering 
A = B = 'i’e and Hi = H e , we find 

= Gfefe + i. 77 PO^) <f G di d x + 2(n PqT) 1 / 2 qQ dl f e . 

(9) 


From Eqs. ©, © with A = B = f e and ©, we 
obtain 


Gfefe = 


npojq-i) 

1 - y/x(q — i) 

(?-*) 


+ 77 Po? 


1 2 


.1 ~y/x{q-i) 
and the mixed Green’s function 


Gd i 


di 


G dlfe = y/^ 1 _%^_ i) Sd ldl , 


( 10 ) 


(ii) 


determined from Eq. © by considering A = g?i , 
B = f e and Hi = H e , with the parameter x = 
(ttPoVi 2) 2 and q = i=i= Here we assume 

the wide band limit denoted by q —> 0 . 

Additionally, for the Hamiltonian of Eq. © we 
have the LDOS p 0 (s) = — ^Im (Gf 0 f 0 ), with 


Gfofo =-^(r)Tr{e b [/ o (r),/t(0)] + } (12) 


and f 0 = Coq■ We notice that Gf 0 f 0 is decoupled 
from the adatoms. Thereby, from Eqs. © and m, 
we take A = B = f 0 and Hi = H 0 in Eq. © and 
we obtain 


~ _ npo(q-i) 

Gfof ° 1 + \fx(q — i)' 


(13) 


Thus the substitution of Eqs. ©, (HU), and (in in 
Eq. ©, leads to 


TOO 

% 


1 + (1 ~ Qb )fl m (^di,di) + 2qbT Re(^di,di)> 

(14) 


where f = HA is an effective adatom 1 -tip cou¬ 
pling, % = H+xy 2 represents the transmittance 
through the STM tips when the adatom 1 is absent 
(f = 0),1Z b = l- T b stands for the corresponding 


reflectance and q b = ^ as the Fano pa¬ 

rameter 44, 45, 47]. The current equation for the 
transmittance encodes three distinct Fano regimes 


of interference as follows: i) q b — > oo (x = 0) pro¬ 
vides T (e) = — Tlm(Gd 1 ,d 1 ) = 7 rTpii, where we use 
pn = — ilm (Gd ly di) as the LDOS of the adatom 
1 , which contains resonant states characterized by 
peaks, since the electronic transport in the inter¬ 
ferometer is solely through the adatom 1 [l5|; ii) 
q b = 0 {x = 1 ) leads to T (e) = 1 — 7 rFpn exhibit¬ 
ing Fano dips as a result of the suppression caused 
by pn over the first term representing the maxi¬ 
mum amplitude 76 = 1. In such a case, when pn 
shows a maximum, T (e) presents a corresponding 
minimum as expected and the tunneling between 
the STM tips becomes the dominant process in 
the system. It reveals the depletion of charge in 
the LDOS of the STM tips detected by the trans¬ 
mittance, once this charge accumulates within the 
adatom 1 , in particular, around the resonant states 
of pn observed in situation (i) for q b —> oo (x = 0); 
iii) corresponds to q b « 0.35 (x = 0.5), which is the 
intermediate case wherein asymmetric Fano line- 
shapes appear with peaks and Fano dips coexisting 
in the same profile, thus making explicit a competi¬ 
tion between V and Vi 2 on an equal footing. These 
features of Fano interference will be addressed in 
SecO 

2.3. System Green’s functions 

By applying the EOM on 

Gdid! = ~j0 (r) Tr{p e [di (r), d\ (0)]+}, (15) 

and changing to the energy domain e. we obtain the 
following relation 

(e- ei- Ti)Gdidi = 1 - tGd 2 ,di - dl (16) 

expressed in terms of the self-energy E = -~{y/x + 
i)r and Green’s functions Gd 2 ,di and G d t d ■ Accord¬ 
ing to the EOM approach we find 


Gd 2 ,di 



and 


tGdi.di ^^d\,dl 

(e — £2 + * 0 +) (e - £2 + * 0 +) ’ 

(17) 

A Gdi,dl 

{e + £2 + * 0 +) (e + £2 + i0 +) 

(18) 


Gd\,di = -2tAKG dl ,d!, (19) 
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which K 
- 


K 


e+ei +T.-K- ' 


with K = 


[ E 2_ E 2_j_ 2iE0+ I_( 0 +)2] ; n as the complex conjugate of 


S and K± = 


[£ 2 -£^+2ie0+-(0+) 2 ] 


Thus substi¬ 


tuting Eqs. (E3- m and m into Eq. m the 
Green’s function of the adatom 1 becomes 


Qd\d\ — 


£ — £\ — X — X] 


MFs 


where 


E MFs = K+ + (2tA) 2 KK 


( 20 ) 


( 21 ) 


accounts for the self-energy due to the MFs con¬ 
nected to the adatom 1. We highlight that the 
self-energy of Eq. m contains the underlying 
mechanism that allows the invariance of the Fano 
lineshape. Such a universal feature is revealed 
as independent on the Fano parameter , which 
we will discuss in detail below. Particularly for 
t = A = we highlight that the expressions for 

K and Emfs found in Ref. [3| are recovered. Such 
a result will be revisited in Sec [3] 


3. Results and discussion 

Below we investigate the features of the system 
Green’s functions by employing the expression for 
the transmittance (Eq. (Thill ). According to Eq. (J3J) , 
this transmittance can be obtained experimentally 
via the conductance G in units of Go for temper¬ 
atures T — > 0. Additionally, we employ values for 
the Fermi energy e, Ej, t, and A in units of the 
Anderson parameter T. 

In Figj2] we consider the Fano regime x = 0 
(i qb —> oo) for the transmittance T of Eq. (Hill as a 
function of the Fermi energy e. This situation cor¬ 
responds to the case where the electron tunneling 
occurs exclusively through the adatom 1, due to the 
strong coupling between it and the STM tips . As 
predicted by the standard Fano’s theory [4J|, the 
transmittance should exhibit a peak around each 
localized state in the adatom probed by the tips: 
see the green line shape of panel (a) for the adatom 
1 here assumed to be decoupled from the adatom 2 
for a sake of simplicity, which leads to the resonance 
centered at £ = E\ = —5 with maximum amplitude 
T = 1 . By keeping this level at such a value and 
employing t = A = 4 combined with £2 = 0, a 
ZBP given by T =1/2 emerges due to the MF ex¬ 
isting in adatom l[5j. Additionally, the most left 


resonance in the same curve corresponds to that at 
E = — 5 found in the green lineshape, in particular 
with renormalized peak position e ss —10 as the af¬ 
termath of the connection y/2A with the adatom 2, 
and with higher amplitude (T > 1/2) in respect to 
that for the ZBP (T = 1/2). Notice that a third 
peak in the vicinity of e « +10 is found charac¬ 
terized by T < 1/2. Thus in the presence of finite 
couplings t. the original peak at e = £1 = —5 in the 
green curve of panel (a) with T = 1 is split into 
those at e ~ —10 and e « +10 both with T < 1 
as the red lineshape points out. In which concerns 
the curves for t < A (t = 2 and A = 4) and t > A 
(t = 4 and A = 2), the transmittance is revealed as 
independent on the strengths t and A. Such a be¬ 
havior attests the situation in which coupled MFs 
are absent in the system within the adatom 2. 



Fermi energy Fermi energy 


Figure 2: (Color online) Transmittance as a function of the 
Fermi energy of the STM tips within the Fano regime x = 0 
{qb —> 00 ): (a) for several cases in the parameters t and 
A with £2 = 0. Particularly in the condition t = A, the 
standard zero-bias peak is found, (b) For t = A the trans¬ 
mittance shares the same lineshape of the symmetric situa¬ 
tions above {£2 = 6T) and below (£2 = — 6r) the MF zero 
mode £ = 0, respectively for positive and negative AFM gate 
potentials. The transmittance then becomes an invariant 
quantity under these conditions. Panels (c) and (d) reveal 
distinct profiles when the system is driven away (t A) 
from the point t = A. 

Panel (b) of Fig(2] depicts the situation in which 
the system is still within the regime t — A, but with 
£2 ^ 0. We point out that such a panel reveals a uni¬ 
versal behavior in the transmittance profile when 
symmetric values for £2 are accounted. This regime 
is characterized by coupled MFs (see Fig.QJd)) 
which result in the suppression of the ZBP and the 
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splitting of the resonances at e ss —10 and e ~ +10 
observed in the red curve of panel (a). For instance, 
the aforementioned universality is verified provid¬ 
ing two identical curves for both values £2 = 6r 
(positive potential) and £2 = — 6r (negative poten¬ 
tial) due to the self-energy Xmfs of Eq. ETT) for the 
MFs, which is dependent on the amplitude K±. No¬ 
tice that K± = [E 2_J t + 2i^o+-(o+) 2 ] within tllis st¬ 
ation, thus implying in -ftT±(£2) = A±(—£2) as well 
as Smfs(£ 2 ) = Xmfs(~£ 2), which then ensure the 
invariance of the transmittance profile at the point 
t = A. For 1 /Awe have K± oc ±£2(t 2 — A 2 ), thus 
allowing a strong dependence on the sign of £2. In 
panels (c) and (d) of the same figure, the transmit¬ 
tance respectively for t < A (t = 2 and A = 4) and 
t > A (t = 4 and A = 2) exhibit distinct behav¬ 
iors as expected when we adopt symmetric values 
£2 = 6r and £2 = — 6r. Therefore in both limits 
t > A and t < A, the influence of the negative po¬ 
tential £2 < 0 on the transmittance is made explicit 
once Emfs(£ 2) y^ X M f s (—£ 2)- 


and A combined with £2 = 0, we can observe 
the crossover from the regime 1 / A (t = 2 with 
A = 4 for the blue lineshape and t = 4 with 
A = 2 in the case of the curve for the orange 
color) towards the point t = A = 4, where we can 
clearly realize in the red curve the emergence of a 
dip with amplitude T = 1/2, analogously to the 
opposite Fano regime of interference (x = 0 and 
qb —> 00) found in Fig. [2(a). In presence of the 
potential £2 y^ 0, the zero-bias dip disappears ac¬ 
cording to the curves with £2 = 6T (positive poten¬ 
tial) and £2 = —6T (negative potential) as found 
in panel (b) of the same figure. As in Fig.[2Kb) we 
also report a universality feature in the transmit¬ 
tance profile, which still arises from the condition 
Xmfs(£ 2 ) = Xmfs(—£ 2) apart from the Fano pa¬ 
rameter as we can notice in Eq. Ell) . For t / A 
coincident curves no longer exist and the universal 
behavior is not verified as pointed out by panels (c) 
and (d), which have the same set of parameters as 
in Fig. [2] 



Fermi energy Fermi energy 

Figure 3: (Color online) Transmittance as a function of the 
Fermi energy of the STM tips within the Fano regime x = 1 
(qb = 0): panels (a)-(d) display the same features of Fig. [2] 
in the opposite regime of interference. 

Fig. [3] holds within the Fano limit x = 1 (qb = 0) 
where the electron tunneling between the STM tips 
is the dominant process in the system, thus result¬ 
ing in Fano antiresonances instead of peaks in the 
transmittance profiles as a function of the Fermi en¬ 
ergy. Panel (a) first displays the case in which the 
adatom 1 is decoupled from the adatom 2 repre¬ 
sented by the green curve characterized by a dip 
at£ = £i = — 5. By using finite values for t 



Figure 4: (Color online) Transmittance as a function of the 
Fermi energy of the STM tips within the Fano regime x = 
0.5 (qb ~ 0.35): panels (a)-(d) display the same features of 
Figs f2l and l3l thus attesting that the point t = A is protected 
against the Fano effect as well as the symmetric change of 
the level in the adatom 2. 

To make explicit that the invariance feature of 
the gate potential for the point t = A = 4 is achiev¬ 
able for any Fano ratio we present in Fig.[4]the 
case x = 0.5 (<#> ~ 0.35) in which both paths V 
and Vi2 of Eq. m compete on an equal footing. 
For this situation, we find intermediate Fano pro¬ 
files where the underlying physics of Figs. [2] and [3] 
is still the same. Therefore based on the results 
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Figure 5: (Color online) Density plots for the transmittance 
as a function of the Fermi energy of the STM tips and the 
potential £2 tuned by the AFM tip in distinct Fano limits: 
(a) x = 0 (qb —»■ 00 ), (b) x = 1 ( q 5 = 0), (c) x = 0.5 
(qb ~ 0.35) for t = A and (d) x = 0.5 (qb ~ 0.35) in the 
situation f / A. The symmetrical panels (a), (b) and (c) 
suggest that a pair of MFs is formed in which only the MF 
A couples to the adatom 1 (see Fig^Jd)). In panel (d), the 
absence of the mirror symmetry under analysis arises from 
the simultaneous coupling of MFs A and B with the adatom 
1. It occurs via the distinct amplitudes ( t + A) and (A — t) 
as Eq. © ensures for t 7 ^ A. 

of Figs 121 [3]and[31 we demonstrate that the invari¬ 
ance of the Fano profile is independent on qb aris¬ 
ing solely from the feature £mfs(£ 2 ) = Smfs(~£ 2), 
due to the connected MFs appearing in the term 
162 ^a'S’b within Eq. 0, in particular for t = A. 
Moreover, the invariance with the potential £2 hr 
the transmittance becomes clearer if we look to its 
density plot spanned by the axes £ (Fermi level) 
and £2. Fig. ©a) is for £1 = —5 and t = A = 4: it 
exhibits the case x = 0 (qb —> 00) for the regime of 
Fano interference, which shows the mirror symme¬ 
try under consideration with respect to the vertical 
axis placed at £2 = 0 (see the vertical dashed lines 
in the same figure). Notice that such a feature also 
manifests itself in panels (b) and (c), respectively in 
the limits x = 1 and x = 0.5. Here the orange color 
designates perfect insulating regions and those con¬ 
ducting are represented by red color. In panel (d) of 
the current figure, this mirror symmetry signature 
is broken just by using t ^ A as expected (A = 2 
and t = 4). 

The invariant Fano profiles for the transmittance 
found in panels (b) of the Figs. [5] [3] and [I] as well 
as those (a), (b) and (c) for Fig. [5] are due to the 


symmetric swap of the adatom 2 level around the 
MF zero mode which reveals the formation of a pair 
of MFs within such an adatom when t = A. It cou¬ 
ples the MF A to the adatom 1 with amplitude \/2A 
as Eq. © ensures. However in the regime t yf A, 
the pair of MFs still exists within the adatom 2, 
but with the MFs A and B hybridized distinctly 
with the adatom 1 via the strengths (t + A) and 
(A — t), respectively. This feature is then probed 
by the symmetric tuning of the potential £2 around 
the MF zero mode, which yields the panels (c)-(d) 
of the Figs© [3] and 31 in addition to the density 
plot in panel (d) of Fig. [5] Therefore, the aforemen¬ 
tioned mechanism ruling the transmittance profiles 
via the couplings of the MFs with the adatom 1 is 
encoded by the self-energy Emfs of Eq. (12T1) . 

4. Conclusions 

In summary, we have explored theoretically a 
Fano interferometer composed by STM and AFM 
tips over superconducting adatoms, in which the 
pair of MFs under the latter, elucidates the gate 
invariance feature of the interferometer, due to the 
Majorana nature arising from the adatom 2. Par¬ 
ticularly for the situation where only one MF hy¬ 
bridizes with the adatom coupled to STM tips, 
the aforementioned invariance consists of a uni¬ 
versal behavior within the transmittance when the 
AFM tip tunes symmetrically the energy level of 
its adatom around the MF zero mode. Such a uni¬ 
versality is constituted by a common Fano profile 
in transmittance as a function of the Fermi level of 
the STM tips for two symmetric values for the AFM 
tip potential. In the case of two MFs connected to 
the adatom beneath the STM tips, we verify that 
such a universality is broken. Hence, despite the ex¬ 
perimental challenging of the proposal, we expect 
that in the near future such an interferometer can 
be developed. 

Acknowledgments 

This work was supported by the Brazilian 
agencies CNPq, CAPES, 2014/14143 - 0 and 
2015/23539 — 8 Sao Paulo Research Foundation 
(FAPESP). 

Appendix A. Transmittance derivation 

Here we use the interaction picture to calculate 
T{e) from Eq.©. It ensures that a state |4>„) from 
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the spectrum of the Hamiltonian given by T~L e + 
Ho — H^-o admits the following time-dependency 


|$ n )= e i/-oo«tun(r)dr |^ n ) 

i r° _ 

~( 1 --/ 'Htun('r)dr) |^„), (A.l) 

hJ-oo 

where H = ^ and |^ n ) is an eigenstate of H v= q. 
Thus the current ^Ttip-i for the STM tip 1 can be 
obtained by performing the expected mean value of 
the current operator Itip-i = Ztip-i (t = 0), which 
reads 


J7tip-1 — (^ro|-^tip-1 |^n) 

i r° 

= / [2tip-i,Htun(r)]dr|® n ) 

+ OtfO, (A.2) 

where we have regarded (din I Ttip-i |^ n ) = 0 and 
by considering the thermal average on the latter 
equation, which gives 

i 

J^tip-i — I Tr{^ = o[Itip-ij ’Htun('r)]}dr, 

(A.3) 

where £>^=0 is the density matrix of the system de¬ 
scribed by the Hamiltonian H v =o- By applying the 
equation-of-motion on Itip-i, we show that 

-7-tip-1 = jT [c ^ ^ C^feClfc ; —0] 

k 

= (“7f ^j V ^ 2 {( C lk d i~ d \ C ek) 

+ (clfcdi - 4 c ofc)} 

^12 ^ ] ( C og C eq — C tq C oq )j 

99 

(A.4) 

which, in combination with Eq. (IA.3I) . leads to 

e f +oa 

Jtip-i = --Ayulm / dr{V2HJ-(-r) 

^ J — OO 

+ 2H 12 A4(-t)}, (A.5) 



= -T d (- T )^{0<P=o[fld 1: J2 c lk( T )c o k{r)}} 

k 

(A.6) 


and 


M(-t) = -^0{-T)Tr{g v= o[flfe,^2cl k {T)c 0 k(T)}} 

k 

(A.7) 


are retarded Green’s functions. 

In order to find a closed expression for the current 
Jiip-i, we should evaluate the integrals in the time 
coordinate r of Eq. (IA.5I) . which result in 


r*+oo 


0 — (3E n p — fiErr 


(e - e 


) 


cItti—t) = z 1 e' ■ 

V ^ E n -E m +i0+ 

mn 

X {ttnl/tdil*™) 

X(^IE4^I^) ( a - 8 ) 


and 


r*+oo 


cLtM{—t) 


Z_1 E 

mn 


( e ~ 0 E n ^ e -(3E m j 

E n - E m + i0+ 


X ^ n \ftfe\*m) 


where we have used Z as the parti¬ 
tion function of H v = ol'hm) = E m |T m ), 
A(t) = e^ Wv=oT Ae~^' hL ' p=oT for an arbitrary 
time-dependent operator A (r). To eliminate the 
matrix element (^ m | c' ek c 0 k |\&n) in Eqs. (IA.8I) and 

(IA.9D . we calculate (^ml Efe c L c °fc^v=o] l^n), 
which gives 


(^ml El c \k c ° k l^«) 

k 


-\/2V(* m \dlf 0 \* n ) 

{E n - E m ) 

2H 12 (^ m \ fjfo\^n) 

( E n - E m ) 

(A.10) 


By performing the substitutions of Eqs. 4MD, 
(1A.9I) with (1A. 101) in Eq. (1A.5I) . we enclose the re- 


where 
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suit into the function labeled by \mn to show that 


jtip-i = \ e ^ z 1 X 

mn 

x 8 (E n - E m ) 


Xn 


(e & En — e P Em ) 

rp _ zt 1 

JZj n - C/ ra 


= -jnA^Y,[ Z ~ le ~ PEn5 ^ E n ~ E m)) 

mn 

X Xnmj (A.11) 

where we have defined 


Xnm = (V2Vf (tf„| /]di |$ m ) <tf m | d\f 0 |tf n ) 

+ 2\/2V ( 2 V 12 ) <tf n | fUl l^m) (^ m | flfo l^n) 

+ (2Fl 2 ) 2 (* n | fife | *m) (*m\ flfo |*n> ■ (A.12) 

In this calculation we have used 


{*n\fldl\* m ) (*m\flfo\*n) 

= (* n \flfe\*m) (^m\d\f 0 \^ n ), 


with 


(e P En — e @ Em ) 

En E m 


—/3e~P En 


(A.13) 


in the limit E n —> E m . The property ['H e ,'H 0 \ = 0 
ensures the partitions E n = E^+E° and Z = Z e Z a 
for the Hamiltonians T-L e and "H 0 , respectively in the 
brackets of Eq. (1A.11I) , thus leading to 


Z~\-^8{E n - E rn ) = J de (~^f) 

x (e -^ +e -^ )(e ^K + e-^) 

x6(s + E e n -E e m )6(e + E°-E° m ). (A.14) 

Therefore, we substitute Eqs. (IA.12I) and (IA. 141) 
in Eq. (IA.11I) to calculate G = = 0). 

The comparison of such a result with Eq. © al¬ 
lows us to find Eq.©. We also verify that G = 
75 ^i7tip-2(V = 0) for the STM tip 2, which is agree¬ 
ment with Refs, jlll. Il2| that show for the case with 
symmetric couplings V absence of Andreev cur¬ 
rents, i.e., j^JTtip-iXv? = 0) J7ti p _2 (</3 = 0). 
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